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Abstract — The aim of the work was to develop a methodology for taking into account the presence of local surface heat 
exchanges in rods of finite length, often taking place in studies of the non-stationary phenomenon of thermal conductivity. 
The proposed methodology was oriented to the subsequent creation of a computational algorithm and its implementation on 
a personal computer using the universal DELPHI programming tool. This allowed the authors to complete the initial stage of 
the research works, which will subsequently take into account the internal heat sources in the rods of finite length and 
constant cross-section in the study of non-stationary thermal conductivity, as well as to develop new methods, approaches 
and models associated with them. 
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I. Introduction 

In mechanical engineering, in particular, in instrument engineering, plastic engineering, and also in other areas, 
physicotechnical processes often arise, where, any structural element in the form of a rod of finite length is suddenly exposed 
to thermomechanical loading [1], Thermomechanical load can be in the form of local temperatures, heat fluxes or heat 
exchanges, where, in the case of their action on rods of finite length and constant cross-section, non-stationary thermal 
conductivity can occur in the system, as well as other non-stationary thermoelastic processes [2], 

In practice, such phenomena are often encountered at the start of gas-generating, nuclear, hydrogen power plants, rocket and 
hydrogen engines, as well as internal combustion engines. In transient, unsteady heat conduction processes in the load- 
bearing elements of these power plants or engines, a complex non-stationary thermally stressed deformed state arises. 
Herewith, in order to study these processes, many load-bearing elements of the above-mentioned structures can be taken as 
rods of finite length with a constant cross-section along the entire length. For an adequate description of a non-stationary heat 
conduction process arising in the rod under the action of dissimilar types of heat sources, taking into account the presence of 
local thermal insulation, it is necessary to use the classical laws of energy conservation [3], because, the application of 
energy conservation laws describing such complex non-stationary processes of thermal conductivity in rods of finite length 
that are exposed to dissimilar kinds of heat sources allows us to take into account natural boundary conditions. As a result, 
the results obtained will have high accuracy [4]. This, in turn, will contribute to a correct evaluation of the thermally steady- 
state behavior of the bearing elements [5], In connection with this, therefore, the development of special effective methods 
for investigating the non-stationary processes of thermal conductivity of load-bearing elements in the form of rods of finite 
length and constant cross-section is an actual problem of the applied non-stationary theory of thermoelasticity. At the same 
time, the development of a complex of applied DELPHI programs that allow studying the classes of non-stationary heat 
conduction processes for the above systems under the influence of dissimilar types of heat sources, taking into account the 
presence of local thermal insulation, is of independent scientific interest. 

II. The Development of the Methodology Results and Discussion 
2.1 Algorithm for the construction of local approximating splines of a function 

Based on the known technique [6] consider a horizontal discrete element of a rod of length l [cm] and a constant cross- 
sectional area F [cm 2 ]. Coefficient of thermal conductivity of the rod material is k x [W/cm°C]. We direct the horizontal axis 




International Journal of Engineering Research & Science (IJOER) 


ISSN, [2395-6992] 


[Vol-5, Issue-10, October- 2019] 


X, which coincides with the axis of the discrete element of the rod. The calculated scheme of the discrete element of the rod 
under consideration is shown in Fig. 1. 

Suppose that the law of temperature distribution along the length of the discrete element is T=T(x), where x is unknown. In 
the local coordinate system ( T,x ), 0 < x < /we introduce the following notation T(x — 0) = 7}; T (x = = Tj; T(x = 

l—Tk. These notations are shown in Fig. 2. 



Element of the Rod Local Coordinate System (T,X) 


We denote the node with coordinate .r=0by/, .r=//2 by j and x =/ by k. n the notation adopted, we can construct the following 
approximating quadratic function 

T{x) = ax 2 + bx + c, 0 < x < l 


a, b, c - are some constant values that are not yet known. In addition, also using the adopted notation, we have 
T{x = 0) = a- 0 + b- 0 + c = 7} 

T(x — /) = a ■ l 2 + b ■ l + c = T k 

Solving this system, we determine the values of the constants a, b, c. 
c — Ti, 


( 1 ) 

( 2 ) 

(3) 

(4) 


al 2 + bl = T k - Ti) 

From the last two equations we ha veal 2 = T k — Tj — bl ; al 2 + 2 bl — ATj — 47’ i ;or7’ fc — Tj — bl + 2 bl — ATj — 47’ i ;orb/ = 
47’ — 4Tj — T k + T ; . Hence we define the value of b 

4Tj-3Tj-T k 


b = 

l 

Substituting (5) into the second equation of system (4), we obtain 
al 2 =T k -Ti + 3 T t +T k - 47} = 2 T k - 47} + 27}; 

From this we find the value a 

2T k -ATj+2Ti 

a — - 5 - 

i 2 

Further, substituting the found values of the constant parameters a, b and c in (1) we obtain 


(5) 


T(x) 


_ 27 ^- 47 -;+ 27-^ 2 + 4Tj-3T t -T k 


-x + Ti, 0 < x < / 


(6) 


(7) 
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We simplify equation (7) according to the scheme(... )7’ i + (...)7J + (...)7j. 

= (2^)7-, + (^)r, + (^)r t ,o < * <, 

Here we introduce the following notation 

.. , ■>. 2x 2 —3 lx+l 2 * t s \ 4lx—4x 2 -. s -v 2 x 2 —lx 

NtW = - ~ 2 -; Nj(x) = ——; N k (x) = ——, 

Then (8) can be rewritten in the following form 

T(x) = N t O) ■ T t + Nj (x) ■ Tj + N k (pc) ■ T k , 0 < x < l 


(8) 

(9) 

( 10 ) 


It should be noted that the values of the node temperatures7), 7}, T fc are still unknown; the functions A/ ; (x), /V ; (x), N k (x)are 
approximate quadratic spline functions. Now we study their properties 

Ni(pc = 0) = 1; Nj (x = 0) = 0; N k (x = 0) = 0 3 

Nt (* = £) = 0; Nj (x = l ~) = 1; N k (x = '-) = 0 1 (11) 

N t (x = 0 = 0; Nj (x = 0 = 0; N k (x = 0 = 1 J 


From (10) it is also possible to determine the expression for the temperature gradient within the length of the discrete element 
under consideration. 


— = ■ T t + ■ 7} + ■ T k = ■ T t + ■ 7) + ^ ■ T k , 0 < x < l 

dx dx 1 dx J dx k l 2 1 l 2 J l 2 k 


In addition, from (9) it can also be determined that 


Ni(x) + Nj (x) + N k (x) 


2x 2 -3lx+l 2 +4lx-4x 2 +2x 2 -lx 

T 2 


2.2 Development of methods for accounting the presence of local surface heat transfer 


( 12 ) 


03) 


Now consider one discrete element, through the lateral surface of which there is a heat exchange with its surrounding 
medium. In addition, through the cross-sectional area of the ends of the discrete element of the rod, heat exchange with the 
surrounding medium also takes place. To begin with, consider the case where the heat transfer coefficient is the same 

everywhere |, the ambient temperature is T ainb [° C], It is required to construct a resolving system of ordinary 

differential equations, taking into account the simultaneous presence of local surface heat exchanges. The calculation scheme 
for the problem under consideration is shown in Fig. 3. 


^ocl * ,l l 


T oc 2> ,l 2 



I 


I 





^oc3> ^3 


X 




1 


T'amb 1 T'crnib 2 1amb 3 f amb 


h t = h 2 — h 3 = h 

FIGURE 3: The Design Scheme of the Problem with Surface Heat Transfer 
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For the discrete element of the rod under consideration, taking into account the presence of local heat exchanges, we write 
the expression for the functional of the total thermal energy 



J , l(T-TjdS + \^{T-TjdS+ J hT-TjdS 

S(x= 0) Z S ls Z S(x=l) Z 


(14) 


Where V is the volume of the discrete element of the rod; S(x=0) is the cross-sectional area of the left end of the discrete 
element; S / s is the area of the lateral surface of the discrete element of the rod; S(x=0) is the cross-sectional area of the right 
end of the discrete element; Jis the coefficient of thermal conductivity of the material of the rod, and X is the heat 

capacity of the material of the rod, whose dimension is A 


In equation (14), the first term 



d\z - physically means the amount of annulled heat quantity in the volume 


of the discrete element of the rod; the second term | [ 


X—~TdV is 4 1C change in the amount of heat per unit time; the 
dt 


other three terms mean the amount of heat that arises from convective heat transfer. Here, it should also be noted that the 
dimension of n is Yl[W ■ °C]- i.e. the value of the heat work performed. In equation (14), the first term - physically means the 
amount of annulled heat quantity in the volume of the discrete element of the rod; The second term is the change in the 
amount of heat per unit time; The other three terms mean the amount of heat that arises from convective heat transfer. Here, 
it should also be noted that the dimension of P is P [W • °C] - i.e. The value of the heat work performed. Now let's proceed to 
integrate all the terms of expression (14). Consider the first term of this functional 



dx. Substituting here the expression (12) we have 


n, 



8N i (x) T | SNj{x) t | dN k (x) 
dx dx 1 dx 


l 

= ~t\ [(4x - 3/)T + (41 - Sx)Tj + (4x-l)T k f dx = 
21 0 

= [lT 2 — 1 6TT + ITT, - 1 6T T. + 1 6T : 2 + IT . 2 \ 

\ 1 1 J 1 K J K J K / 


Here the sum of the coefficients will be zero (7-16+2-16+16+7)=0. 
Further we integrate the second integral of the functional (14) 


c dT r oi 

TF = \ A—TdV = FM— Tdx = 
2 J At J At 


ST 


dt 


dt 


0 

:FAj 


dr 


dT i dT, 

Al.(x)-i + AC(x)-^ + AF(x)-^ 
dt dt dt 


■ N,(x) ■ T, + (x)-Tj+N k (x) ■ T k ]ix = 


dT dT dT 

N t 2 (x)-Tr-L + N,(x) • Nj(x )• Tj ■ -4 + N,(x)■ N k (x)■ T k ■ —4 + 
dt dt dt 


= N, (x) ■ Nj (x) ■T l ~ + N ] 2 (x) ■T r S ^ + Nj (x) ■ N k (x)-T k -^ + 


dt 


dt 


st st st 

+ N, (x) • N k (x)-T r —^ + Nj (x) • N k (x) • T } ■ + N k \x) ■ T k ■ —4- 

dt dt dt 


dx 


(15) 


(16) 


In this equation, integrating each integral (of nine) separately and substituting them into (16), we obtain the integrated form 
ofH 2 ; 
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n 2 = fai 


1 dr l dr i dr 1 dr j 

j _ 1 \ T — !_ — T — 1 l T _— 

15 ' dt 15 j dt 30 A ' dt 15 ' dt 


8 „ dT l dT l dT k 1 dT k 2 dT k 

—Tj — — H- T k — - T - 1 - T - 1 - T k - 

15 1 dt 15 dt 30 dt 15 J dt 15 dt 


Now, in the expression (14) we calculate the third integral 


n 3 = 


S(x= 0) 


j f l(T-TjdS=™(Ti-Tj = ^{ 


Flu 


= ^-{T 2 -2T i T oc +T oc 2 ) 


Similarly, we can calculate the fifth integral in expression (14) 

n.= f T-TjdS = ^{T k -Tj = ^ 


—\T 2 -2T k T oc + Tj) 


S(x=l ) 


In the expressions (18) ... (19) the sum of the coefficients will be equal to zero (l-2+l)=0. 
The computation of the fourth integral in expression (14) will not be easy 


n 4 = J hT-rJds = ! ^\(r-rjdx = ^J[dv,.(x) -t + n/x)- r. + N k (x)-T k )-rJdx= 


Ph 


Ph 


s, 

Ph r 
~2 


l 

\[(N i (x)-T i +N j (x)-T j +N k (x)-T k )-T oc l(N i (x)-T i +N j (x)-T j +N k (x)-T k )-T oc }lx = 


Ph r 


(17) 


(18) 


(19) 


( 20 ) 


J k 2 W • T i + 2N M) • Nj(x) ■T r T j + 2 IV, (X) • IV, (x) • T ■ T t - 21V,. (x) • T • T oc + 

~ 0 

+ N 2 (x) ■ T 2 + 2Nj (x) • N k (x) • r. • T k -2N j {x)-T r T oc + N 2 {x)-T 2 -2 N k (x)-T k -T oc + T oc 2 \lx 

Where P is the perimeter of the cross-section of the discrete element of the rod. Now calculate each integral in expression 
( 20 ) 


K l = l \N 2 (x)-T 2 dx = ^T 2 


( 21 ) 


K 2 = j 21V,. (x) • Nj (x) • r ■ Tjdx = ^ TT j 


15 


K 3 = j 2N, (x) • N k (x) ■T r T k dx = ~- 7)7) 


( 22 ) 


(23) 


K, 


2 r 


= | - 21V,. (x) • 7). • T oc dx = - J (2x 2 - 3/x + / 2 )7 ]T oc dx 


2 ,2x 3 3/x 2 , , 

■ — (-+ /x) 

Z 2 3 2 


o 


?7 3 4_9 + 6 / 

TT =——(-- —-)TT = -—TT 

i oc 1 2 v /- ' i oc n IOC 

L 6 3 


(24) 


t: 


8/ r 2 


\n 2 (x)-T 2 dx- 

0 1 J 

ZV 6 = j 21V. (x) • IV, (x) • Tj ■ T k dx = TT k 


(25) 


(26) 
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K 7 = | (- 


I / 

j (-2 Nj (x) • Tj ■ T oc )dx = --\ (4 lx - 4x 2 )TjT oc dx = 


2 2 4x 3 2/ 3 6-4 4/ 

- - y~) T i T oc =-~r (—) r , T 0C = - y T i T oc 


(27) 


z 1 ' 1 z 

^8=1 IV t 2 (x) • 7) ~dx = |(2x 2 - lx)T k 2 dx = —^\ (4x 4 -4/x 3 + l 2 x 2 )T k 2 dx= 

n In l n 


1 4.T 5 4 /V 

= -f(-zx +-) 

Z 4 5 3 


T 2 = — (12-15 + 5)7) 2 =— T 2 
15 15 


(28) 


a : 0 


2 r 


■ \{-2N k (x)-T k .T oc )dx = -^j(2x 2 - lx)T k T oc dx 


0 

2 2jc 3 Zjc 2 

I 2 ( 3 2 } 

/ 


2/ 4-3 / 

T k T oc =--r(—)T k T oc =--T k T oc 


K 


to 


Jr^=r. 2 .-ad' 


oc lo 


IT: 


6 


o 


(29) 


(30) 


Now substituting (21)... (30) into (20) we obtain the integrated form n 4 


n 4 = 


Phi 


'2 2 2 1 1 8 2 2 4 2 2 1 2l 

_ t- 4 - tt - TT—TT h-T _ h- TT - TT -t - T~—TT +T 2 

15 , i5 ,J 15'* 3 ' oc 15 2 15 2 * 3 2 oc 15 k 3 k oc oc * 


(31) 


It should be noted here that the sum of the coefficients before the temperatures will be zero 

2 2 1182421 

15 + 15 _ T5 _ 3 + 15 + 15 _ 3 + T5 _ 3 + 1_0 

Substituting (15), (17)... (19) and (31) into (14), we find the integrated form of the total thermal energy functional: 
Fk, 


n = - 


61 


-{it. -1 67)7} + 2T,T k - 167)7) +1 6T 2 + lT k )+ 


+ FM\ 


I of 1 dT. 1 ^ 8T 1 dT 

—I - 1 —T - T, - 1 —I -h 

15 dt 15 2 dt 30 dt 15 8t 


+ 1 t 8j ± + 1 t 8 JL-± T d A + l T d Ir + ^ T d A. 

15 1 dt 15 * dt 30 ' dt 15 2 dt 15 * dt . 
+ — (t 2 -2TT + 7’ 2 )+ — (r, 2 -27)7’ +T 2 )+ 

r\ \ l l OC OC ) rt \ k k OC OC J 


(32) 


Phi 


2 2 118 2 4 2 1 l 

—T 2 + —TT- — TT --TT +—T 2 + —TT-Ztt + — T ~ — Tt t + t 2 


15 


15 


15 


15 


15 


15 


Further, minimizing the functional 77 with respect to the node values 7), 7) and 7) we obtain a resolving system of first-order 
ordinary differential equations with the corresponding initial conditions and taking into account existing local heat 
exchanges. 

dn Fk r , _ , _ „ . I 1 8T t 1 dT 1 dT k 


_ = °. ^ (147) - 167) + 27)) + Fzl Z ( ) + 


Fh Phi f 4 2 

+ T (27)-2T oc ) +-(-7, + -7). 


15 dt 15 dt 30 dt 
= 0 ; 


1 1 \ 

15 Tk 3 Toc ) 
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377 Fk x , . „ / 1 37) 8 dT 1 dT k \ 

——- = 0; -> —f (-167; + 327] - 16T t ) + F/l? ——^ ^ + 

dT. 61 v 1 ' kJ V15 dt 15 dt 15 dt 1 


Phi / 2 16 2 4 

“2“ (15^ + 15 7 '' + 15 _ 3 Tc 


,) = 0 ; 


377 Fk x , . „ / 1 dT t 1 37; 2 dT k 

w t = 0: “ -6f ( 2r ‘ - 167 i +14,, ‘) + ^' (- 30 af + Is af + Is-ar 1 + 

Fh Phi ( 1 2 4 1 \ 

+ — (27k - 2 T 0C ) +— T t + —T + —T k --T oc = 0. 

2 k oc 2 V 15 1 15 7 15 3 oc 7 

Or, after simplification from the last system, we get: 


40 ky 


3J\ 3Tj_ _ 1ST*. I35k^ 15ft 2Pft , 

3t 3t * 3t U / 2 FAj ‘ VfA Al+ 


-V,-+ 


T J+ [s± + J*- 
\Xl 2 2 fA 


dTj 0 dTj d T k I Ph 40k x \ T t80k x 8Ph | ^ / Ph 80k x \ ^ 

st st st va ‘ Uf f!/ 7 ' va k 


I rji _ 15hT 0C SPhT 0C 

\ lb —-1- 

A1 2 f A 

10 PhT„ r 


a 


-l a JA + a 21 + 2^1+ ( 5fc * Ph | r , ( ZL_l£k]7. + / 35k x 15 h 2 Ph \ r _ 5P/ir 0 , 

i3f 3f 3f 1 1 - 11 ; 1 1 1 J 7 Uf Ai fa; fe fa 


(33) 


Setting the initial conditions 

7) lt=o = a i! Tj\ = a 2 ; 7’ fc | t=0 = a 3 

one can find a solution of the system (33). 

The obtained results are oriented for the computational algorithm created by the authors, which was realized in a personal 
computer in the form of Delphi programs [7], With the help of the Delphi-7 programming developed on the object-oriented 
programming language, we solved the non-stationary heat conduction problem for a rod of finite length under the influence 
of local heat transfer along the lateral surface [8], 

III. Conclusion 


Based on the energy conservation laws, a technique is developed that simultaneously takes into account the presence of local 
surface heat exchanges for the non-stationary heat conduction problem occurring in a rod of finite length and the constancy 
of its cross-sectional area. Applying the quadratic approximation function of the form from the energy conservation law by 
the method of minimizing it through the node temperature values, we obtain solving systems of first-order linear ordinary 
equations that take into account the natural boundary conditions. Due to this, the results obtained have a high degree of 
accuracy. The developed methodology and the corresponding computational algorithm allowed to realize calculations in a 
personal computer on the object-oriented Delphi-7 programming language, and thus to solve the non-stationary heat 
conduction problem for a rod of finite length under the influence of local heat exchange along the lateral surface. 
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